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Abstract. We present in this paper a study of a sample of 14 nearby clusters of galaxies (0.06 < z < 0.1) 
observed with the ROSAT/PSPC. We only select clusters with low galactic absorption (nH< 6 x lO^^cm"^) in 
order to trace the hot X-ray emitting intra-cluster medium (ICM) out to large radii. We convert the X-ray surface 
brightness profiles of the clusters into emission measure profiles scaled to the classical scaling relations based on 
the spherical collapse model. We sort the clusters into different temperature bins and stack the scaled emission 
measure (ScEM) profiles of clusters belonging to the same bin together. This method enhances the statistics of 
the profiles - especially in outer regions. The stacked profiles allow us to observe a signal out to radii r > r2oo. 
In the center (r < 0.4r2oo) we find that the ScEM profiles deviate from predicted scaling laws. Hotter clusters 
have systematically a higher ScEM than cooler clusters. This result is in very good agreement with current studies 
on the Lx — T relation and the entropy-temperature relation {S (x yO-^S) found recently. At radii r > 0.4r2oo 
we find that the ScEM profiles agree well within the error bars, suggesting self-similarity. Fitting beta-models to 
the overall ScEM profiles we find for the different sub samples Vc = 0.15 — 0.18r2oo and (5 = 0.8, which is higher 
than the canonical value of /3 = 2/3 often found. The beta-model is generally a better representation for hotter 
than for cooler clusters. We see indications for continuous steepening of the profiles with increasing radius: at 
radii r > 0.8r2oo the profiles are systematically below the beta-model curve with /3 = 0.8. 

We discuss our results with respect to the observed X-ray luminosity Lx — T relation, the gas mass Mgas — T 
relation and the total mass M — T relation. We also address implications on the origin of the observed S ~ T 
relation. Furthermore we discuss the observed steepness of the X-ray profiles, which falls off more rapidly than 
predicted from the NFW-profile for cold dark matter halos. 
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1. Introduction 

In the hierarchical scenario of structure formation, clus- 
ters of galaxies are the largest and youngest objects in 
the universe which have virialized. This makes them ideal 
targets for cosmological studies. In order to use clusters 
as probes for cosmology one has to be sure that internal 
processes did not have time to erase the primordial infor- 
mation from which they have formed and that gravity was 
the only driving force responsible for the formation of the 
cluster. 

X-ray observations of clusters of galaxies allow to 
study cluster physics using the hot intra-cluster medium 
(ICM) , which is spread out all over the cluster (see Sarazin 
119861 for an introduction). In the past years studies have 
shown that the central parts of clusters show the pres- 
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ence of non-gravitational physical processes such as cool- 
ing flows or cores (e.g. Fabian '1994'2UUi5' ), interactions 
with AGN's (e.g. Boehringer et al. il993i ; Belsole et al. 
120011 McNamara et al. 12005(1 or sources, which change 
the ICM temperature from the virial temperature in the 
central parts (see Voit l2()()^ for an overview). Therefore, in 
order to study cosmology with clusters it is wiser to con- 
centrate more in the outer parts of clusters, where internal 
cluster processes are yet largely inactive due to much lower 
densities, which require longer timescales. Unfortunately 
lower density goes together with lower emission: the ICM 
emits energy mainly via thermal bremsstrahlung, which 
is proportional to the density squared. Therefore exter- 
nal cluster regions show less X-ray emission with respect 
to the center, where the density is much higher. This 
lower emission translates in lower statistics for available 
X-ray observations. There are several ways to overcome 
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the problem of low statistics in outer regions. One pos- 
sibility, which is often used is to fit the central parts of 
clusters with models, typically the beta-model or one of 
its more sophisticated derivatives, and to extrapolate into 
outer and undetected regions. The disadvantage of such a 
method is clear: the extrapolation depends on the validity 
of the model in external cluster regions. In order to over- 
come this problem we choose here to go a different way. In 
order to enhance statistics we employ a stacking method, 
in which we add different cluster profiles together. Before 
we add the profiles together we have to make sure that 
they are similar. For this we apply general cluster scaling 
laws which have in the past shown to work for observed 
X-ray cluster profiles (e.g. Mohr et al. 119991 Vikhlinin et 
al. tl999: Neumann & Arnaud 'MJT; Arnaud et al. 

1^2,) and translate the X -ray profiles into emission mea- 
sure profiles, which represent the gas density squares in- 
tegrated along the line-of-sight. 

Since we want to trace clusters out to large radii 
we need X-ray observations with large field-of-view 
and low background. Therefore we choose to look at 
ROSAT/PSPC (TriimperUnnSl observations of a sample 
of nearby clusters. 

The paper is organized as follows. After the introduc- 
tion we describe briefly the sample selection in Sec. 2. In 
Sec. 3 we review the scaling relations which we use for the 
cluster profiles and the stacking. Sec. 4 describes the data 
treatment. In Sec. 5 we present our results, which we dis- 
cuss in Sec. 6. We finish with the conclusion in Sec. 7. 

We use the following cosmological parameters through 
out the paper: ri,,„o = 0.3, I^a = 0.7, Hq = 70 km/s/Mpc. 
All error bars are 1 a. 



2. The sample 

2.1. General selection criteria 

In order to trace the cluster profiles up to large radii we 
choose to use ROSAT/PSPC data which offer a large field- 
of-view of 2 degrees. These data have further the advan- 
tage of showing low background contamination. In order 
to assure that the cluster emission is well traced out to 
the virial radius with a large observed region only having 
background emission, which is important for background 
subtraction, we restrict ourselves to the redshift interval 
0.06 < z < 0.1. Furthermore, to avoid high absorption 
of X-ray emission, which is especially critical in low sur- 
face brightness regions, such as the outskirts of clusters, 
we choose clusters with a hydrogen column density of 
nH < 6 X lO^^cm^^. We only consider clusters with a 
fitted temperature estimate (preferentially obtained with 
ASCA or Beppo-SAX). Tab^ shows the clusters in our 
sample with some of their physical properties. In total we 
selected 14 clusters for our analysis in the temperature 
range 1.7 <kT< 8.5 keV. 



2.2. Notes on cluster temperature selection 

When there exist several kT measurements for one clus- 
ter we take the latest published measurement. In most 
of the cases the references come from Ikebe et al. (2002). 
Comparing the results obtained by these authors with oth- 
ers, specifically with Markevitch et al. 119981 we find gen- 
erally good agreement within the error bars. 

There exist several kT measurements of A2670, which 
range from 3.7 (Novicki et al. IMM White, Jones & 
Forman TW) up to 5.6 keV (Sanderson et al. HUU^ . 
There are several intermediate measurements which agree 
within the error bars at around 4.5 keV (Horner et al. 
1999: Finoguenov, David & Ponman^^), so that we de- 
cided finally to include this cluster in our sample at this 
temperature. There were two clusters, which in principle 
fulfill our criteria, but which we discarded for further anal- 
ysis: the first cluster is Abell 2734. Ikebe et al. lMT^ fitted 
kT=5± 0.4 keV for this cluster but Reiprich & Bohringer 
"2002 estimated kT=3.9 ± 0.6. These different results are 
incompatible and we therefore exclude this cluster for our 
study. The second cluster, which we skipped is RXJ2344. 
The reason: Ikebe et al. 120021 give large uncertainties of A 
kT=4 keV, which we consider too high for our analysis. 

3. Scaling relations 

3.1. Virial radius R and r2oo 

Assuming spherical symmetry the mass of the virialized 
part of a cluster is related to its virial radius by 



M/R^ = 



(1) 



(e.g. Kaiser IT986fl Defining Ac as the density contrast 
with Ac = Pvir/Pcriu ~ Pcrit = 3iJo/(87rG) is the critical 
density of the imiverse and E'^(z) ~ Pcrit{z) / pcrit{z = 0) 
we can write: 

MIR' ^ A,p,„tiz) oc Ac P"-;'^'^\ = A,E'{z) (2) 

Pcr'i.t(Z = 0) 

Applying the virial theorem to the ICM gives M oc 
RT, where T is the ICM temperature. This yields together 
with eq[21the mass temperature relation 



M cx 



y3/2 



^cE{z) 

and the radius temperature relation 

Roc ^ , , 
^cE{z) 

with 

E'^{z) = noii + zf + n^ 

in our choice of cosmological parameters and 
Ac = IStt^ + 82{n{z) - 1) - 39{Q{z) ~ \f 



(3) 

(4) 

(5) 
(6) 
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Table 1. General cluster properties. For the references: I : Ikebe et al. I2(JU2I G: Gomez et al. 119971 D : De Grandi & 
MolendiEOna H: Horner et al.^UM 



name 


nH in 




kT 


Ref. 


7^200 


7^200 


expo, time 




10^" cm-^ 




in keV 




in Mpc 


in arcmin 


in s 


A578 


4.41 


0.0870 


1 7+^-'" 

'-■ ' -0.3 


G 


1.12 


11.8 


4519 


A644 


5.14 


0.0704 


6.5 ± 0.3 


I 


2.23 


28.2 


10246 


A1651 


1.71 


0.0845 


6.2 ± 0.5 


I 


2.15 


23.1 


7429 


A 1750 


2.49 


n 0852 


4.5 ± 0.3 




1.83 


19.6 


13146 


A 1775 


1.00 


0.0716 


3.7 ± 0.3 


I 


1.68 


20.8 


8669 


A 1795 


1.20 


0625 


6.2 ± 0.3 


I 


2.19 


30.8 


62076 


A2029 


3.07 


0.0773 


7.9 ± 0.4 


I 


2.45 


28.4 


15693 


A2142 


4.05 


0.0909 


8.5 ±0.5 


I 


2.50 


25.2 


17215 


A2244 


2.07 


0.0968 


5.8 ±0.6 


I 


2.06 


19.6 


2963 


A2255 


2.51 


0.0806 


r q+0.4 


I 


2.10 


23.6 


14534 


A2597 


2.50 


0.0852 


4.2 ±0.5 


I 


1.77 


18.9 


7164 


A2670 


2.92 


0.0762 


4.5 ±0.2 


H 


1.85 


21.7 


17679 


A3112 


4.90 


0.0750 


4.7 ±0.3 


I 


1.89 


22.5 


7598 


A3921 


2.80 


0.0936 


4.9 ±0.4 


I 


1.90 


18.7 


11997 



with 



n{z) 



E^{z) 



(7) 



using the usual spherical top hat model (Eke et al. 
119961 Bryan & Norman 1998 ) . In the case of a critical 
density universe, is constant with IStt^. 

For the normalization of the x-axis of the surface 
brightness profiles we use r200j the radius in which the 
mean over-density of the cluster is 200 times the critical 
density of the universe. Numerical simulations have shown 
that within this density contrast clusters can be assumed 
to be in equihbrium (see Cole & Lacev ITM?)jl . At larger 
radii in-fall becomes important and the assumption of viri- 
alization is not anymore valid. To calculate r2oo we use its 
definition: 



M200 = 



47r 



-200pcrit''200 



(8) 



We are interested in the r2oo-T relationship at the 
mean redshift of our sample, z=0.08. With the choice 
of our cosmological parameters we find Pcrit = 146 
/kpc*^. To calculate r2oo we need a relationship between 
M and r.To establish this relationship we can use the hy- 
drostatic equation using the X-ray emitting intra-cluster 
medium (ICM): 



1 dP 
Pa d?- 



d$ 
dr 



GMlr) 



(9) 



Pg is the gas density and P is the pressure. We can 
use the equation of state of a perfect gas for the ICM: 
P — nkT. Assuming that the ICM is isothermal we only 
need the gas density distribution of the gas to calculate the 
total mass at a given radius r. Referring to the beta-model, 
which we will introduce in Sec l5.3l and more specifically in 
eqI221we can rewrite eqElas: 



M{r) = 



3fc/3 Tr^ 



/i = 0.59 is the mean molecular weight. We find for /S 
values of 0.80 (see Sec l5.3|) . Since we are interested in the 
outer parts of clusters, where r >> Tc, we can neglect Tc 
for getting r2oo- 

Therefore we can approximate at large radii: 



M(r 



3fc0.8Tr 



Gijmp 

Inserting eg II II into eqlHl gives: 



kT 

7-200 = 0.86\/— Mpc 



(11) 



(12) 



at z = 0.08, which is our r2oo-T relation at this red- 
shift. We will discuss this and the resulting M — T relation 
in detail in the discussion. 

Now, practically, to scale the profiles to r2oo we do the 
following: we calculate for each cluster the angular diam- 
eter distance and translate the surface brightness profile 
into a profile with Mpc scale using our adopted cosmolog- 
ical parameters. We then apply the relationship shown in 
eq^for each cluster to place each cluster at z=0.08.- We 
multiply the radius in Mpc scale by {yJ /S.c{z = 0.08)£^(z = 
0.08))/(-\/Ac(z)£^(z)). Please note: the redshift interval is 
very small so that this correction is not very important. 
Now, that we "virtually red-shifted" the cluster to z=0.08 
we can apply ea ll2l to calculate its r2oo and divide the Mpc 
scale of the surface brightness profile by this value. The 
surface brightness profile is now in units of r2oo- In Tab^ 
we display the physical parameters which are necessary to 
calculate r2oo for each cluster. 

3.2. Emission measure profiles and self-similarity 

The emission measure is defined to be the gas density 
squared integrated along the line-of-sight: 



(10) EM{r) = 



nl{l)dl = 



^ n^{x)xdx 



(13) 
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Assuming self-similarity, the mass-temperature and 



radius-temperature relation and supposing that M, 



gas 



M we find: 



EM{r) oc oc \/T(VAI^(z))3 



(14) 



however, EM{r) is not measured directly. What can 
be observed is the X-ray surface brightness of ICM (see 
below), which is 



S{r)=K{T,z) 



EM{r) 
47r(l + z)i 



(15) 



A(T, z) is in our case the emissivity in the chosen 
ROSAT band (0.5-2.0 keV) taking into account interstellar 
absorption and the PSPC spectral response. A(T, z) can 
be calculated by simulating a red-shifted cluster spectrum 
with a given plasma temperature and folding it through 
the spectral response of the ROSAT/PSPC. In our case 
we use XSPEC for the spectral simulations. One can thus 
translate the observed surface brightness profile into the 
emission measure EM{r) profile. Self-similarity of clusters 
implies that the scaled emission measure profiles should be 
identical for all clusters: 



ScEM{r/r2oo) 



47r(l + z)45(r/r2oo) 



K{T,z)^T{^,E{z)f 



(16) 



ScEAI{r /r2oo) can be calculated directly from the ob- 
served surface brightness distribution. We apply eg II 61 for 
all the cluster profiles using the best fit temperature dis- 
played in TabQ 

4. Data treatment 

4.1. Surface brightness profile 

For each cluster we create an image in the 0.5-2.0 keV 
range and calculate the corresponding exposure maps tak- 
ing into account wobbling of the satellite and the geome- 
try of the support structure of the PSPC. We create from 
these images radial surface brightness profiles with a bin 
size of each annulus of 15 arcsec. Obvious point sources 
or clearly identified substructures are cut out. All profiles 
have 220 bins (or 55 arcmin) with the exception of A1775: 
this cluster was observed off-axis and we have to restrict 
ourselves to 190 bins in this case. 

We use the software package EXSAS for these tasks 

4.2. Background subtraction 

Since we are interested in the profiles of the outskirts of 
the clusters where source emission is low, background sub- 
traction is a crucial point in our analysis. 

To determine the background for each cluster we cal- 
culate the mean emission Bi in the radial range 1.5 < 
T/r2oo < 2.5. In this region we are confident to not any- 
more detect cluster emission. In order to check whether the 
background is flat we also calculate the background in the 



region 1.5 < r/r2oo < 2 (S2). If |Bi - 2 x a/y/no.\ < B2, 
where a is the standard deviation of B2 and no. is the 
number of radial bins which fulfill the criterion 1.5 < 
f 1^200 < 2, we consider the background emission as flat. 
There are some clusters (A644, A2029, and A1795) for 
which our proflle does not go up to 2 x r2oo- In this case 
we consider for Bi the region in which r/r2oo > 1- and 
for B2 1.0 < r/r2oo < 1-5. All clusters fulfill the crite- 
rion of I Bi — 2 X a I ^prio.\ < B2 with the exception of 
A3112 and A2670. In these particular cases we determine 
the background Bi in the interval 1.0 < r/r2oo < 1-5. 
For all clusters we use Bi as background estimate in the 
following, which we subtract from the original extracted 
surface brightness profile. The error of the background es- 
timates is typically 10% of the statistical error in the sur- 
face brightness bins close to r2oo- Adding this background 
uncertainty quadratically to the statisical error would en- 
hance the error by roughly 1%, which is negligible. In the 
following we therefore do not take into account the error 
of the background. 

We applied several methods to determine the back- 
ground level and checked the robustness of our results 
with different background subtractions. All background 
subtractions we performed yielded the same results. 

4.3. From the surface brightness profiles to the scaled 
and stacked emission measure (ScEM) profiles 

To obtain the ScEM profile for each cluster we use the 
background subtracted surface brightness profiles and ap- 
ply ea ll6l to the y-axis. The scaling to r2oo is explained in 
SecO 

4.4. Stacking and sub sample selection 

We are interested in the outer parts of clusters, where 
statistics is generally sparse. In order to improve the statis- 
tics we opt for a stacking method, in which we add dif- 
ferent emission measure profiles together. We also want 
to study the temperature dependence of the scaled emis- 
sion measure profile. Therefore we divide our sample up 
in sub samples which represent different cluster temper- 
atures. FigEl shows the temperature distribution in our 
sample and the selection of the 6 sub samples in different 
temperature intervals. The sub sample composition is also 
shown in Tab [3 

FigQlshows the scaled emission measure profiles for all 
clusters before stacking. 

In order to stack the scaled profiles together all the 
points of the different ScEM profiles must be located at 
the same relative position with respect to r2oo , the radius 
to which we scale the profiles. For this we interpolate all 
scaled emission measure profiles linearly with a step- width 
of 0.03 X r2oo- This gives for each profile 33 points up to the 
virial radius. This is about half of the points we have for 
each original emission measure profile. With the reduction 
of the data points we make sure that each original data 
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Table 2. The membership list of clusters of the different sub samples. The first column shows the sub sample number, 
the second column the temperature interval for each sub sample and the third column shows the average temperature 
of each sub sample. 



no. 


kT 


kT 


name 
















min-max 


[keV] 


([keV]) 
















[keV] 


















1 


1.7-4.9 


4.0 


A578 (1.7) 


A1775 (3.7) 


A2597 (4.2) 


A1750 (4.5) 


A2670 (4.5) 


A3112 (4.7) 


A3921 (4.9) 


2 


3.7-4.9 


4.4 




A1775 (3.7) 


A2597 (4.2) 


A1750 (4.5) 


A2670 (4.5) 


A3112 (4.7) 


A3921 (4.9) 


3 


4.2-4.9 


4.6 






A2597 (4.2) 


A1750 (4.5) 


A2670 (4.5) 


A3112 (4.7) 


A3921 (4.9) 


4 


5.8-8.5 


6.7 


A2244 (5.8) 


A2255 (5.9) 


A1795 (6.2) 


A1651 (6.2) 


A644 (6.5) 


A2029 (7.9) 


A2142 (8.5) 


5 


5.8-6.5 


6.1 


A2244 (5.8) 


A2255 (5.9) 


A1795 (6.2) 


A1651 (6.2) 


A644 (6.5) 






6 


7.9-8.5 


8.2 












A2029 (7.9) 


A2142 (8.5) 



A 

P- 
g 

p. 
1 
1 

..— t 


mi 



10-5 




r/r200 



Fig. 1. The scaled emission measure profiles of all clusters 
before stacking. The error bars here are based on pure 
photon statistics. 



6 

5 h 



0) 



X) 
3 



4 
3 
2 
1 h 



23456789 
cluster temperature in keV 

Fig. 2. The temperature distribution of the clusters (see 
also Tab^ in our sample and the selection in temperature 
of sub sample 1 to 6. 



point is only used once for the linear interpolation and 
that the resulting data points are thus uncorrelated. We 



calculate the mean scaled emission measure for each sub 
sample: 



: ScEMkir/r2oo) 



J + 1 



(17) 



k is the index for the sub sample and ranges from 1 
to 6 (see also Tab|2l. J ^^^^ n is the number of the first 
and of the last cluster in the sub sample. The number of 
the cluster raises with temperature, so the coldest cluster, 
A578 has number 1 and the hottest cluster, A2142 has 
number 14 (see also TabEJ. We also calculate ScEM{r) 
for all 14 clusters together. 

4.4.1. The uncertainty of the scaled emission measure 
profiles 

There are different possibilities to calculate the uncer- 
tainty of the scaled and stacked emission measure profiles. 
As a first attempt one can add quadratically the different 
statistical errors of the individual profiles based on photon 
numbers. 



AstatScEMk{r/r2ao) = 



J2-^^{As,,,ScEM{^,^W 
n — j + 1 



(18) 



A.BtatScEM{r/r2oo) is the mean error of the stacked 
profile taking into account the internal error of each pro- 
file. However, in this error calculation differences linked 
to effects such as morphology (substructure, ellipticity) , 
redshift or temperature uncertainties are not taken into 
account. Another approach is to calculate the dispersion 
for the stacked profiles at a given radius and to determine 
the mean error from this. The mean error per point for 
the stacked profiles is in this case: 



AScEM{r/r2oo) 



cr{r/r2oa) 



(19) 



where (T(r/r2oo) is the rms at r/r2oo averaged over all 
profiles in the sub sample. 

If only uncertainties due to photon statistics were 
important AstatScEM{i) and AScEM{i) should be of 
the same order. Else, AScEM(i) has higher values than 
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AstatScEM (i). Therefore, to be conservative we choose in 
the following AScEM{i) for the error bars. As we will see, 
in the central parts there is a systematic trend of higher 
ScEM with higher temperature. This increases the disper- 
sion of the scaled emission measure profiles for each sub 
sample and is also, of course, dependent on the tempera- 
ture interval of the sub sample. The smaller the temper- 
ature interval, the smaller should be the apparent disper- 
sion. This effect is not of statistical origin but is systemat- 
ics. This results in a likely overestimate of the calculated 
error bars in the central parts. 

5. Results 

5.1. The different stacked emission measure profiles 

Fig|2| shows the scaled emission measure profile averaged 
over all clusters as well as the mean profile of sub samples 
1 and 4 (see also TablSJ, which each represent half of all 
clusters. Up to about 0.4r2oo we see a systematic trend 
of higher ScEM for higher kT clusters, which signifies a 
deviation from self-similarity and implies that apart from 
the spherical collapse based on gravity, additional clusters 
physics plays a role in these regions. At larger radii the 
different ScEM profiles agree well within the error bars. 
In Fig2|we show the stacked emission measure profiles of 
the smaller sub samples. Also here, in the central parts 
the same trend is visible. The higher the average cluster 
temperature the higher the central values for the ScEM 
profile. The profiles start to agree within the error bars 
at different radii. This effect might be at least partially 
caused by the low number of clusters in the sub samples. 
Sub Sample 6 does not allow a proper assessment of the 
statistical uncertainties since it only contains 2 clusters. 
Additionally, one of the clusters in this sub sample is Abell 
2142, which is known to contain an important amount of 
substructure (Markevitch et al. I2UUU|I . 

We show in Fig|Slthe scaled emission measure at differ- 
ent radii as a function of average sub sample temperature 
(see TabEJ. Displayed is at the same time the rms of the 
best fit cluster temperatures divided by the square root 
of the number of clusters in each sub sample (error bars 
in x-direction). In Fig|5lwe see again up to about 0.4r2oo 
a systematic trend of higher ScEM for hotter clusters, 
which disappears, when approaching 0.4r2oo- We will dis- 
cuss the differences of the ScEM profiles in the central 
parts in SecEl 

5.2. Quantification of differences 

In order to quantify the differences between the different 
profiles we apply a x^-test in different radial intervals with 
ri < r < r2 , in which x^ is defined in the following way: 

2 _ (S^EW.in) - ScEMiin))^ 
AScEMi (r, ) + AScEMi (r, ) 

AScEMi is defined in eqE| We only concentrate here 
on sub sample 1 and 4 since they are the only ones in which 



T3 
0) 

"3 
O 

m 




10-9 ^ 



<kT> = 6.7keV (4) 
<kT> = 5.4keV 
(all clusters) 
<kT> - 4.0keV (1) 



10-1 



10° 



r/r200 



Fig. 3. The stacked scaled emission measure profile of sub 
sample 1 and 4 and for all clusters together. 



Table 3. x test results between sub sample 1 and 4 in 
different radial regimes. The first two columns give the 
inner and the outer boundary of the radial interval taken 
for the test. 



int. radius 


ext. radius 


x' 


d.o.f. 


X' / d.o.f. 


in r2oo 


in r-200 








0.1 


0.2 


16.9 


2 


8.45 


0.1 


0.3 


35.4 


6 


5.9 


0.1 


0.4 


43.3 


9 


4.81 


0.1 


0.7 


50.0 


19 


2.61 


0.1 


1.0 


67,5 


29 


2.33 


0.2 


0.3 


18.5 


3 


6.17 


0.2 


0.4 


26.4 


6 


4.4 


0.3 


1.0 


32.0 


22 


1.46 


0.4 


1.0 


24.1 


19 


1.27 


0.6 


1.2 


22.8 


19 


1.20 


0.7 


1.2 


20.8 


16 


1.30 


0.8 


1.2 


19.8 


13 


1.52 



we have enough clusters to achieve sufficiently good statis- 
tics for this test. The corresponding results are shown in 
Tabini The results vary from 8.5 for the central parts 
down to 1.2 in the external parts of the clusters. The lower 
the value of the higher is the probability that the 2 pro- 
files have the same shape and normalization. The lower 
values for the external regions are explained by the larger 
error bars and the apparent missing trend that hotter pro- 
files show a higher scaled emission measure. It is a recon- 
firmation that the profiles do not coincide in the central 
parts. 



5.3. Beta-model fits 

The beta-model is defined the following way see (Cavaliere 
& Fusco-Femiano J976J: 



^(r) = 5o(l + ^) 



-3/3+0.5 



(21) 
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Fig. 4. The mean scaled emission measure profile and the 
mean profiles for sub sample 1, 2 and 3 (top panel) and sub 
sample 4, 5, and 6 (bottom panel). Please note: sub sample 
6 comprises only 2 clusters and the shown error based on 
the dispersion of the profiles has to be interpreted very 
cautiously. 



Fig. 5. The scaled emission measure as a function of tem- 
perature for different radii. The corresponding radii are 
shown in units of r2oo right to the figure. The error bars 
in the x-direction show the statistical uncertainty of the 
cluster temperature in each sub sample (see text). The 
horizontal lines show the mean ScEM averaged over all 
the sub samples at the selected radii. The dotted line 
shows the relationship ScEM oc T and the dashed line 
ScEM (X T^-^. Please note: sub sample 6 comprises only 
2 clusters and the shown error based on the dispersion of 
the profiles has to be interpreted cautiously. 



S{r) is the surface brightness. 5*0 is the central surface 
brightness, Tc the core radius and (3 a slope parameter. 
The beta-model is such a successful model since it allows 
to deproject easily the surface brightness of the ICM which 
emits via thermal bremsstrahlung into a gas density pro- 
file: 

n(r) = no(l + 4)"'^/' (22) 
r 

' c 

Emission measure profiles have the same shape as sur- 
face brightness profiles. Therefore we can use the beta- 
model also for our stacked emission measure profiles. To 
further quantify the differences of the various scaled emis- 
sion measure profiles we fit beta-models in different ra- 
dial intervals. Again, for statistical reasons we only con- 
centrate here on sub sample 1 and 4. Tab0] resumes our 
results. We fit in 4 different radial intervals. We do not 
look at regions < 0.1 x r2oo since central parts of clusters 
are often affected by so-called cold cores, which raise the 
central surface brightness of some clusters and which sub- 



sequently enhance the dispersion of the profiles in these 
regions. 

As can be seen in Tab0]for the fit in the different re- 
gions the parameters P and Vc change less for sub sample 
4 (the hotter clusters) than for sub sample 1 (the cooler 
clusters). This indicates that the beta- model is a better 
description for hotter than for cooler clusters. In FigElwe 
show the ScEM profiles with logarithmic binning and the 
best fit beta-model profiles for 0.1 < r/r2oo < 1.2 . We 
see that the beta-model is too shallow for sample 1 in the 
center and at radii r > 0.8r2oo- For sample 4 we see that 
the beta-model is a good description in the center, how- 
ever, also in this sample we see indications that the profile 
falls off more rapidly in the outskirts than predicted with 
our best fit overall beta- model. In Fig[7|wc zoom into the 
external parts of the ScEM profile and show the results of 
the different sub samples with logarithmic binning. There 
is an obvious trend that the data points lie below the 
beta-model curves at radii r > 0.8r2oo- 
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Table 4. Results of the beta-model fitting. Error bars are 
ItT-level. The first column shows the sub sample number. 
The second shows the interval in units of r2oo in which the 
beta- model is fitted, is also in units of r2oo- 



sa. 









5*0 X 
10"'' 


d.o.f. 


Area. 


1 
1 
1 
1 


0.1-1.2 
0.1-1.0 
0.1-0.7 
0.3-1.2 


o.i8g:i? 
o.o89mt 
1.4^:^4 


o.8o!i:r, 
o.78jJ:?i 
o.64[S:i 


2.1 
2.3 
6.6 
2.5 


33 
26 
16 
26 


1.42 
1.66 
0.48 
0.72 


4 
4 
4 
4 


0.1-1.2 
0.1-1.0 
0.1-0.7 
0.3-1.2 


o.i5t!:{| 

0.136g:l|g 

o.i3g:il 

U.^Oo,22 


0.77g:?^ 
0.89^! 


5.8 
6.5 
6.9 
1.8 


33 
26 
16 
26 


0.68 
0.61 
0.30 
0.88 



In Fig|B| we see that the fitted beta-models for sub 
sample 1 and 4 are identical at radii r > 0.6r2oo- This can 
be easily explained by the fact that we fit identical f3's 
for the 2 samples. The differences found for the core radii 
do not play anymore an important role in the outskirts of 
clusters. 

The found Xred values are systematically lower for sub 
sample 4 than for sub sample 1. This again indicates that 
the beta-model is a better description for hotter than for 
cooler clusters. 

For the two samples we find a mean /3 of about (3 = 0.8 
and — 0.15r2oo for sub sample 4 and r^T = 0.19r2oo for 
sub sample 1. The mean value for /3 which we find is with 
around 0.8 higher than the canonical value of 2/3 as often 
found by fitting individual clusters without any stacking. 

If we perform beta- model fits with f3 — 2/3 fixed, we 
obtain values of Vc of 0.09r2oo for sub sample 1 (the cooler 
clusters) and 0.039r2oo for the hotter ones. The corre- 
sponding reduced values are 1.9 (sub sample 1) and 
1.6 (sub sample 4). These values are significantly higher 
than the reduced values for our best fits leaving /3 as a 
free parameter. 

5.4. Power law fits 

Beside the beta-model fits we also perform power law fits 
in different radial intervals for sub sample 1 and 4. The 
power law wc fit has the following form: 

S{r) = N2oor-" (23) 

A^200 is the normalization of the power law at r2oo and 
a is the power law index. Our fit results are displayed in 
Tab[Sl We fit power laws with increasing a as we go to 
larger and larger radii, which shows a continuous steepen- 
ing of the profiles. In the outskirts of the cluster a beta- 
model with P — 0.8, as found in the previous section gives 
a = 3.8. The power law fits we find for radii > 0.7r2oo are 
systematically higher than this. This indicates again that 
the ScEM profiles show a steeper slope in the external 
parts of the cluster than predicted by a beta-model with 
P = 0.8. a = 3.8 is not compatible with any of our fit 
results in the external parts of the cluster. 




]^Q-10 I I I I I I I I lU 

10-1 loo 

r/r200 



Fig. 6. The profiles of sub sample 1 and 4 with logarith- 
mic binning with different models superposed. The curve 
corresponding to a = 4.3 is our lower limit for the slope in 
the cluster outskirts (see text). This slope matches at the 
same time a density profile which follows the NFW-profile 
at these radii (see text). We normalize this a = 4.3 profile 
to match the data points at r = 0.7r2oo- 

10-"^ F 1 1 1 1 \ ^ 




]^Q-10 I \ \ \ lJ \ Ll 

loo 

r/r200 

Fig. 7. Zoom into the outskirts of the ScEM-profiles. The 
curve showing a power law with slope a — 7 .22 is based 
on the fit in the outer parts of all ScEM profiles stacked 
together (see also TablSJ. The symbols and error bars of 
the different sub samples are shifted with respect to each 
other only for visibility. The points of the different profiles 
have all the same location in x-dircction. 

The power law fits in the regions 0.7 < r/r2oo < 1-2 
suffer from the fact that we a) have relatively large error 
bars per bin and b) do not have a large number of bins. 
This explains our low Xred values and the large error bars 
for the fit parameters. Since sub sample 1 has such a high 
value for a we also show the fit result for sub sample 
2, which contains all clusters of sub sample 1 with the 
exception of A578, the coolest one. We see a dramatic 
difference of the fitted a in sub sample 1 and 2. There is 
a good agreement between the fit results of sub sample 
2 and sub sample 4 within the error bars. However, sub 
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Table 5. Results of the power law fittmg. Error bars are 
ItT-level. The first column shows the sub sample number. 
The second column shows the interval in units of r2oo in 
which the power law is fitted. 





fmax 


i V200 ^ 




d.o.f. 


Area. 


1 


n 1 1 

U.l-l.Z 




^- '^2.73 


Q/i 
o4 




1 


0.1-1.0 


^•^'-'5. 58 


2 7fin lA 


27 


2.71 


1 


0.1-0.7 


7.87?-?^ 




17 


0.65 


1 


0.3-1.2 


4.i5t:II 


Q or3.52 

'5-'J'-*3.18 


27 


1.87 


1 


0.7-1.2 


o.5og:^? 


11.4^*7? 


14 


0.47 


2 


0.7-1.2 




6.99i"2^^ 


14 


0.51 


4 


0.1-1.2 




3-032i98 


34 


1.58 


4 


0.1-1.0 


c 71=6.28 


s.ooi:^! 


27 


1.46 


4 


0.1-0.7 


6.79^:^1 


2.89i:if 


17 


1.03 


4 


0.3-1.2 


n r74.45 


3.79t;| 


27 


0.80 


4 


0.7-1.2 


oq3.77 
^•002.05 


5.734 47 


14 


0.71 


all 


0.7-1.2 




7 ooW.Ul 
' -^^5.62 


14 


0.438 



sample 1 and 4 give, as sub sample 1 and 2 results not 
compatible within the error bars. We also show the power 
law profiles fitted in the outer parts of the clusters in FigEl 
and Fig[7| All power law profiles show a steeper slope than 
the beta-model. 

The error bars on the ScEM profiles are very large in 
the external parts and it is therefore impossible to give an 
exact value for a and to estimate with certainty whether 
the different profiles have the same shape, e.g. whether 
clusters are 100% self-similar in the external parts or not. 
We see very steep profiles in the outskirts. Because of 
the steepness of the profiles, slightly wrong temperature 
measurements, which shift the profiles along the x-axis 
can have dramatic effects on the power law fit results. 
Currently, we are only able to give a lower limit on a at 
radii r > 0.7r2oo which is a < 4.3. This lower limit is 
based on the fact that it is the lowest value, which is com- 
patible with the fit results (based on the fit of sub sample 
2). 

6. Discussion 

6.1. The beta-model parameter /3 as function of 
temperature 

We fit systematically higher /3-values to the overall cluster 
profiles than previous studies, which find generally /3 = 
2/3. Additionally, previous work has shown an increase 
of P with temperature (e.g. Castillo-Morales & Schindler 
Sanderson et al. lWiSjl .^ 
Previous studies extracted and fitted individual X-ray 
profiles. The disadvantage of looking at profiles without 
averaging them before fitting is the limited statistics at 
the cluster outskirts. However, this part is essential to 
measure the right slope of the profiles. In Tab01we show 
our results of beta-model fitting in different regions of the 

^ A dependence of /3 on the temperature changes the slope 
of the M — T relation, which we will discuss in Sec l6.4l 
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clusters. We can see that the fitted /3's are systematically 
smaller in the central parts in comparison to fits going out 
to r > r2oo- This is especially true for sub sample 1, which 
comprises the coolest clusters in our sample. (This trend 
is in agreement with the results of Bartelmann Steinmetz 
11996! who fit beta-models to the ICM of simulated clus- 
ters) . We therefore suggest that the apparent (3 — T cor- 
relation is linked to the fact that the beta-model is not 
a good overall description of the clusters, especially for 
cold and small systems and that fitting only central parts 
of clusters gives systematically lower values for f3 than 
overall beta- model fits. Since cooler clusters offer gener- 
ally less statistics than hotter systems, due to lower lumi- 
nosity and emission measure {EAI oc VT), they are gen- 
erally detected out to smaller radii with respect to r2oo 
than their hotter counterparts. It is therefore not surpris- 
ing that studies based on fitting individual clusters find 
a trend of increasing (3 with cluster temperature. We can 
resume for clusters at the same redshift and with the same 
exposure time: the hotter the cluster, the larger the radius 
of detected cluster emission with respect to ^200, and the 
higher the corresponding fitted /3- value. 

6.2. The Lj;-T Relation 

The central part of clusters is the region where most en- 
ergy is radiated. Therefore, the central parts of clusters 
are responsible for the majority of the X-ray luminosity 
Lx of clusters. Using a beta-model with rc = 0.15r2oo 
and (3 = 0.8, based on our best fit results, we calculate 
that 90% of Lx comes from the region r < 0.4r2oo- We 
see deviations from the predicted self-similarity in regions 
smaller than 0.4r2oo- At larger radii the ScEM profiles of 
sub sample 1 and 4 agree well within the error bars, which 
suggests self-similarity in the outer cluster parts. 

Previous work (Edge & Stewart lEMl Arnaud & 
EvrardHnSI Markevitch lT^ Allen & Fabian ITOUSjl has 
shown that the observed Lx-T relation deviates from pre- 
dictions based on simple self-similarity. Instead of 

Lx^T\^MlMT)/R') (24) 



Lx oc (25) 

with a certain error was found. It is logically to assume 
that the deviations from self-similarity in the central re- 
gions of clusters and the observed Lx — T relation are 
linked together. If we assume that clusters are self-similar 
for r > 0.4r2oo (which implies Lx in these parts), 

we can calculate the needed relationship Lx oc for 
r < 0.4r2oo to obtain the observed overall Lx oc rela- 
tion: 

0.9T"^ + O.IT^ = (26) 

this gives x ~ 3.1 — 3.5 for temperatures which vary, 
as in our sample between 2 and 8 keV. If our hypothesis is 
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correct, we should observe a Lx ~ T relation within a pro- 
jected radius r < 0.4r2oo which scales as Lx oc j"3 i ■■3-5^ 
Assuming that R (x T^^^, and that the differences origi- 
nate entirely from the Mgas — T relation at r < 0.4r2oo 
with 

Mg,,(r< 0.4r2oo)«r^ (27) 

we find using eq[21 y = 2 — 2.2 instead of y = 3/2, 
which is predicted from self-similarity. 

Our sample ranges from roughly 2 keV up to 8 keV 
clusters, y = 2—2.2 implies that a 8 keV cluster has 2-2. 6x 
more gas mass than predicted from self-similarity with 
respect to a 2 keV cluster((8/2)2-3/2 = 2, (8/2)2-2-3/2 ^ 
2.6 ). However, this is only valid for r < 0.4r2oo- The 
fraction of gas mass in this region is about 6-7% of the 
total gas mass calculated up to r2oo- The increase of Mgas 
is therefore very modest for the entire cluster. 

We can check whether this predicted Mgas — T relation 
in the central region matches our data. ScEM should be 
independent on the temperature in the self-similar case 
and ScEM cx T^ - -^-^ to explain the observed Lx - T 
relation. FigEl shows ScEM oc T and ScEM cx T^-^ for 
different radii normalized to 6 keV. We see indeed, that 
the curves match quite well the data points in the cluster 
center, which reinforces the idea, that it is the observed 
deviation from self-similarity in the central parts of clus- 
ters, which explains the Lx — T relation. In another rep- 
resentation we show in Fig|Sl the ScEM divided by the 
mean temperature for each sub sample. The different pro- 
files should therefore coincide if they scale with T and the 
resulting gas mass with Algas oc instead of M cx T^/^^ 
which is predicted from simple scaling laws. As one can see 
at r < 0.4r2oo the profiles coincide rather nicely, however, 
at larger radii the profiles diverge. 

The idea that the deviation of the gas mass from 
self-similarity expectations is responsible for the observed 
Lx — T relation is not new, we presented this idea already 
in another paper (Neumann & Arnaud 1200 1|. However, 
we see for the first time, that this break of self-similarity 
is only present in the central parts. 

6.3. Entropy in the cluster centers 

It is now common to write the entropy of the ICM in terms 
of S* cx T/rie^^. If clusters were perfectly self-similar the 
ICM density should be independent on temperature and 
thus 5 cx T should be observed. However, S cx T°'^^ is 
found (Ponman et al 2.003 Pratt & Arnaud 1^0051 Pratt & 
Arnaud f2005|l . We observe in the central parts of the clus- 
ters in our sample that ScEM cx T instead of a ScEM 
independent of temperature. The deviation of S from self- 
similarity and the observed ScEM-T relation represent the 
same fact: n cx VT. This result is in agreement with our 
previous work (Neumann & Arnaud l2001|l . where we stud- 
ied a sample of clusters in the z interval 0.04 < z < 0.06. 
The deviation of the entropy from the self-similar predic- 
tion is commonly an argument for the existence of non- 
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Fig. 8. The different ScEM profiles divided by the mean 
cluster temperature for each sub sample. This corresponds 
to Mgas oc instead of M cx T^^^ predicted from simple 
scaling laws. 

gravitational heating in clusters even though recently ra- 
diative cooling to explain the feature became an issue (e.g. 
Muanwong et al. l^nU^ . 

We can synthesize the results: the observed Lx — T 
relation of Lx oc T^, S* cx T^-^^ and ScEM cx T are 
caused by the fact that n cx \/T instead of a density which 
is independent on the temperature in the central parts of 
clusters at r < 0.4r2oo- 

Despite this good agreement of the work by differ- 
ent authors, there exist some important discrepancies: 
Ponman et al. 1)2003(1 observe the entropy scaling of S* oc 
yo.es very large radii, which is in clear contradic- 

tion with our results, which suggest self-similarity in the 
outer parts (r > 0.4^200) and thus S oc T. Ponman et 
al. (2003) fit beta- models to the surface brightness pro- 
files to individual clusters in their sample and then add 
together profiles of clusters in different temperature inter- 
vals^. This is, at first sight, similar to our procedure here, 
however the main difference is that we add the profiles 
together before we fit the beta-model and not after fitting. 
Adding the profiles together before fitting has the advan- 
tage of enhancing the statistics of the profiles themselves, 
especially in our case in the outer cluster parts. These re- 
gions generally provide very low statistical quality with 
respect to the center. Only in adding up different profiles 
we have sufficient statistics to trace the X-ray profiles out 
to r2oo and beyond. 

The differences of the S — T relation in the center and 
in the outer cluster parts is an important piece of informa- 
tion for the physical modeling of the ICM. A lot of effort 
has been already put into explaining the observed devia- 
tion of the S — T relation from self-similarity. Cavaliere et 
al. ll999l developed a model in which hierarchical clustering 
and thus in- fall from substructures was taken into account. 
Tozzi & Norman (2001) and Dos Santos & Dore ipHTI^ 

^ See also Sanderson & Ponman I2UU3I 
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looked at the effects of an initial background entropy level 
and shock heating. More recent models based on sim- 
ulations integrate galaxy feedback and radiative cooling 
(Kay et al lMTl Kay 120031 Muanwong et al l^UH^ . Voit & 
Ponman 1)2003^ propose smooth accretion models to ex- 
plain the observed entropy scaling. For a comprehensive 
summary on entropy profiles and their impact on clusters 
as well as general cluster properties see Voit l)2004|l . Very 
recently another approach was presented by taking into 
account additionally AGN and quasar activity (Lapi et 
al. l2UU5|l . which yields good agreement with observations. 

The fact that we only find departures from self- 
similarity in the inner cluster parts suggests that the 
source(s) responsible for the relation S oc r^-^s act(s) only 
or mainly in the center of clusters. One can see that the 
ScEM profiles are parallel in the centers (see Fig|3| and 
FigEI) (see also Neumann & Arnaud f 19991 Pratt Arnaud 
EOOHl- This imphes that 5 oc T'' '^^ and n oc ^/WF not only 
at r = 0.1r2oo as suggested by Ponman et al. I2UU3I but 
all the way up to r = 0.4r2oo- The offset of the parallel 
curves with respect to each other defines the deviation of 
the S ~T relation from the expected simple scaling laws. 
The fact that the ScEM profiles are parallel suggests that 
at each radius the deviation of the S—T relation is strongly 
correlated to the gas density n. It is therefore logical to 
assume that it is a process which depends on n which is 
responsible for the observed S — T relation. A natural pro- 
cess in this respect would be radiative cooling. Radiative 
cooling could also explain why we only see the departures 
from self-similarity in the center. In the outer parts the 
density is not high enough to make this mechanism sufh- 
ciently efficient to be observed. 



6.4. The M' 



200 



T relation 



Making use of the beta-model and taking eq^] and eqlHl 

together, we can calculate the r2oo ~ T and M200 ^ T 
relation at any given redshift: 



'"200(2:) 



and 



M2oo{z) 



0.89 / kT 



;Mpc 



E{z) V keV 



.0 X 1013 ( kT\ 



(28) 



3/2 



(29) 



Please keep in mind that these are the relations as- 
suming that the ICM is isothermal. 

Bryan & Norman l|1998|l calculated a. M — T relation 
of clusters based on Eulerian hydrodynamic simulations. 
Using our choice of cosmological parameters they found 
for an over-density /S.cE'^{z) = 101 at z=0: 



M = 12.4 X 10 



\keY J 



3/2 



Ma 



(30) 



We can calculate our M — T relation for the same over- 
density (multiplying the Af2oo ^ T relation based on the 
beta-model with ^200/101), which gives: 



M = 11.2 X 10 



13 



/ kT 
IkeV 



3/2 



Mr. 



(31) 



We see a difference of 10% with respect to the results 
of Bryan & Norman H1998|l . The corresponding observed 
R — T relation has a 3% lower normalization than the re- 
lation based on simulations. This difference is quite small 
and is less than found in other studies based on X-ray 
data, which give normalisations much lower than the re- 
sults of Bryan & Norman (1998 J (e.g. Horner et al. 119991 
Mohr et al. 119991 Finoguenov et al. 120011 Nevalainen et 
al. lMMH Allen et al lMTTl Arnaud et a\ TM^ . ^ The dif- 
ference between our results and previous studies is the 
fact that we see a steeper slope for the profile in outer 
parts, which translates in a higher value for (3. Using the 
hydrostatic equation, the assumption of isothermality of 
the ICM and the beta-model to calculate the M — T rela- 
tion one finds that the normalization of this relation is di- 
rectly proportional to (3. Our results of a steepening of the 
ScEM profiles is in agreement with a study of Vikhlinin et 
al. H1999|l . These authors fitted beta- models in the outer 
parts of clusters and found systematically higher /3's with 
respect to fits performed taking the cluster centers into 
account. 

We have seen that the slope of the ScEM profile in 
the outer parts is quite likely steeper than the beta-model 
with (3 — 0.80. We can calculate the ICM density profile 
and subsequently the M200 — T relation taking the power 
law model (ea l23|) used to fit the surface brightness profile: 



n{r) oc r"("+i)/2 



(32) 



/ I rri \ 3/2 

M200 = 7.4 X 10^2 X (a + 1)3/2 / _ \ (33) 

As one can easily see, a higher value for a translates 
in a higher normalization of the M — T relation. 

We observe that the ScEM profile drops more rapidly 
than the beta-model with (3 = 0.8 at radii larger than 
0.8r2oo- Power law fits in the outer parts give a power law 
index a > 4.3. Taking a — 4.3 we find a normalization of 
the M — T relation of 12.6 (see ea l33|l . which is in very 
good agreement with the value of 12.4 found by Bryan 
& Norman l|1998|l . We will come back on the issue of the 
slope in the next section. 



6.5. Indications for instabilities at large radii 

Recently, analyzing a sample of nearby clusters observed 
with XMM-Newton Pointecouteau et al. 120051 (see also 



Please note: when the M — T relation is not given for M200, 
we extrapolated this relation assuming a constant density gra- 
dient and isothermality 
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Pointecouteau et al l2l)U4|l have found, that the matter pro- 
files in clusters fit very well the NFW-profile (Navarro, 
Frank & White [TUMI I1997|l . which has the density distri- 
bution 



p(r) = 



Po 



(r/rs)(l + r/rs 



(34) 



The results found by Pointecouteau match results 
based on Chandra data (Arabadjis et al. 120021 f2()()4l Buote 
& Lewis EHOl Lewis et alEOISJ- 

There is general agreement between observations and 
simulations (Dolag et al. 2004) that the concentration pa- 
rameter c = r2QQ/fs is approximately 6, which implies 
that rg = 0.16 — 0.17r2oo- We show in Fig|H|the density 
distribution of the NFW-profile, assuming = 0.16 and 
compare it with a beta-model slope with (3 = 0.8, and a 
slope with r^**-^^, which corresponds to the density profile 
inferred from the best fit power law in the outskirts av- 
eraged over all ScEM profiles (see TabEJ. We also show 



FigEli 



-2.65 



which corresponds to the density fall-off 



of the ICM if 5 = r with a = 4.3, which is our lower 
limit at r2oo (see also Sec l6.4|l . As one can see in this fig- 



ure, the slope of 



is parallel to the NFW density 



profile around r2oo- Therefore it is not unlikely that the 
ICM density falls off more rapidly close to r2oo than pre- 
dicted from the NFW profile for cold dark matter. How 
can we interpret this result? 

At some radius in the outskirts of clusters a strong 
temperature gradient must occur, since we know from 
COBE results (Smoot et a ll991|l that the material sur- 
rounding a cluster has a temperature of roughly 3 K. The 
ICM drops from a few 10^— 10® K down to a few K within a 
short radial range. This large temperature gradient makes 
the ICM convectively unstable and as a result the cold 
material from the outskirts and the hot ICM mix. This 
mixing is enhanced by accretion of cold material from the 
cluster outskirts. This very likely causes that the ICM at 
large radii is not anymore in hydrostatic equilibrium. It is 
very likely that our result reflects this scenario. We only 
see the hot ICM emitting in X-rays, but we do not see 
cooler material or in general gas, which is below a tem- 
perature of about 10^ K. It is therefore likely that the total 
baryon fraction in the outskirts is similar to the one ob- 
served closer to the cluster center. However, at larger and 
larger radii the hot phase becomes less and less important 
and cooler components not emitting in X-rays contribute 
more and more to the total baryon budget. It is currently 
not clear how the different phases mix in detail and how 
effective heat conduction is in this respect. More detailed 
studies on cluster outskirts similar to the study presented 
here are necessary to give quantitative numbers. For this 
it is important to have the temperature information of the 
ICM at radii >> 7'5oo- Unfortunately, this is currently very 
difficult to obtain. One would need a very high through- 
put X-ray observatory with good energy resolution, low 
instrumental background and large field-of-view (FOV) to 
do this kind of study. XMM-Newton and Chandra unfor- 
tunately do not have sufficiently low internal background 



0.01 - 



to 




0.001 



r/r200 

Fig. 9. Different density profiles in arbitrary units nor- 
malized to the same value at r = ?'200- Since the NFW 
profile and the profile representing r^^'^^ have such an 
identical shape we choose a small offset between the two 
curves at r2oo to enhance the visibility of the plot. The 
line r~^'^^ represents the shallowest density profile, which 
is in agreement with our power law fits on the ScEM pro- 
files in the external cluster parts. The curve with r^^-^^ 
represents the density distribution which originates from 
the best power law fit to the ScEM profile averaged over 
all clusters in the outer regions with ScEM{r) oc r~^-^^ 
(see also TabEj). 

at high energies and provide only a limited FOV to fulfill 
these requirements. Maybe more distant clusters which fit 
in the FOV of XMM-Newton and Chandra can be used 
for a stacking method to constrain the temperature dis- 
tribution of the ICM at large radii. 

7. Summary and conclusion 

In the study presented here we looked at a sample of 14 
nearby clusters of galaxies in the redshift range 0.06 < z < 

0. 1 with low galactic absorption for which good tempera- 
ture measurements are available in the literature. For each 
cluster we calculate the emission measure profile scaled to 
r2oo and apply the spherical collapse model (Eke et al. 
ITOMI and Bryan & Norman ITOUS). We divide the clusters 
in sub samples which represent different temperature in- 
tervals. In order to enhance the statistics, especially in the 
cluster outskirts, we stack the profiles of clusters belonging 
to the same sub sample together. Comparing the different 
scaled and stacked emission measure (ScEM) profiles we 
find the following results: 

1. The scaled emission measure profiles are not self- 
similar in the center, however, in a logarithmic repre- 
sentation the different profiles are parallel in the range 
0.1?'2oo < r < 0.4r2oo- The difference of the profiles can 
be approximated with ScEM oc T which implies a gas 
density n oc VT. This is in agreement with entropy 
studies, which find generally S oc T°-^^ (Ponman et 
ali2003i Pratt & Ar naud 1200112001 ■ The fact that the 
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different ScEM profiles are parallel up to r < 0.4r2oo 
suggests that the difference of the S—T relation, which 
defines the offset of the profiles with respect to each 
other is strongly coupled to the ICM density itself. The 
only currently recognized mechanism which links the 
ICM density to the observed S — T relation is radiative 
cooling. 

2. We find that about 90% of the X-ray emission of clus- 
ters comes from the region encompassing r < 0.4r2oo- 
Therefore this central region is primarily responsible 
for the observed X-ray luminosity. The deviation of the 
ScEM profiles from self-similarity in the centre which 
is responsible for the observed S — T relation explains 
also the observed Lx — T relation Lx oc T^ (Edge & 
Stewart P91, Allen & Fabian T?liJg, MarkevitchTOT, 
Arnaud & Evrard ri999|l . which deviates from simple 
expectations based on structure formation and gravity. 

3. At radii r > 0.4r2oo the different cluster profiles 
agree well within the error bars, which suggests self- 
similarity. This implies that global cluster properties 
defined at large radii, such as the Af — T or i? — T re- 
lation follow the theoretical scaling of M c>c T^/^ and 
R oc Vt. 

4. Fitting beta-models to the overall stacked profiles we 
find generally (3 — 0.8, which is higher than the canon- 
ical value of P = 2/3, which is often observed. The 
found higher value of /3 translates into the need of a 
higher normalization of the M—T relation with respect 
to previous studies. Assuming isothermality up to large 
radii we find a normalization of the M — T which is 
only 10% lower than the relation based on numerical 
simulations found by Bryan & Norman H1998fl . 

The beta-model is generally a better representation 
of the hotter than of the cooler clusters. We find rel- 
atively poor agreement between the beta-model and 
the ScEM profiles in the central parts (r < 0.15r2oo) 
of cool clusters. Simple scaling laws predict that the 
emission measure scales with EM cx This effect 
causes that cooler clusters have systematically smaller 
detection radii with respect to r2oo than hotter clus- 
ters. The bad representation of the X-ray profiles of 
the central parts of cool clusters and EM cx \/T are 
quite likely the reasons why previous studies looking 
at individual cluster profiles found a correlation of in- 
creasing f3 with increasing temperature. 

5. The ScEM profiles are systematically below the best- 
fit beta-model at radii r > 0.8r2oo, which implies a 
steepening of the X-ray profiles with radius. Fitting 
power laws to the external parts of the ScEM profiles, 
we find steeper profiles than predicted from the NFW- 
profile, which fits observed cluster density profiles well 
up to rsoQ. We argue that this steepening of the ob- 
served profiles is caused by the fact that at radii close 
to r2oo the ICM is not anymore in hydrostatic equi- 
librium, as predicted from numerical simulations. At 
these large radii the X-ray emitting ICM is mixed very 
likely with cooler material accreted from larger radii. 



This multiphase gas causes that part of the baryons 
are not anymore detected in X-rays. 

Our cluster sample spans over the temperature range 
1.7keV< kT < 8.5keV, however only one cluster has a 
temperature below 2 keV and only one cluster is above 
8 keV. Therefore our dynamic range in temperature is 
limited. We can therefore currently not exclude that the 
behavior of small groups with temperatures well below 
1 keV or above 9 keV is in agreement with our results or 
not. 
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